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ABSTRACT 



Context. Hydrogen peroxide (HOOH) was recently detected toward p Oph A. Subsequent astrochemical modeling that included 
reactions in the gas phase and on the surface of dust grains was able to explain the observed abundance, and highlighted the importance 
of grain chemistry in the formation of HOOH as an intermediate product in water formation. This study also predicted that the 
hydroperoxyl radical HO2, the precursor of HOOH, should be detectable. 
Aims. We aim at detecting the hydroperoxyl radical HOt in p Oph A. 

Methods. We used the IRAM 30m and the APEX telescopes to target the brightest H0 2 lines at about 130 and 260 GHz. 
Results. We detect five lines of HOt (comprising seven individual molecular transitions). The fractional abundance of HO2 is found 
to be about 10~ 10 , a value similar to the abundance of HOOH. This observational result is consistent with the prediction of the above 
mentioned astrochemical model, and thereby validates our current understanding of the water formation on dust grains. 
Conclusions. This detection, anticipated by a sophisticated gas-grain chemical model, demonstrates that models of grain chemistry 
have improved tremendously and that grain surface reactions now form a crucial part of the overall astrochemical network. 
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1. Introduction 

Water is an essential molecule in star-forming regions, in par- 
ticular because it is the main constituent of the icy mantles 
of dust grains, one of the main repositories of oxy gen, and in 
certai n conditions an important gas coolant (e.g. iNisini et"al] 
2010). It is moreover detected in virtually all star- forming en- 
vironments in which it i s searched for (e.g. ICaselli et al.ll20ldl 
Hog erheiide et al.ll201 ll Ivan Dishoeck et aU201 lb . 

Water is believed to form efficiently on the surface of dust 
grains through several different pathways. The relative impor- 
tance of these pathways is not completely understood, however, 
and the reaction barriers involved are still under debate despite 
numerous recent l aboratory experiments (Iop polo et al.l 2010; 
Cuppen et al. 2010). This is because isolating a particular re- 
action in the laboratory is extremely difficult, and the absolute 
measure of the reaction rates and barriers on the grains is almost 
impossible. 

In this sense, astronomical observations of the precursors of 
water can be of great help to constrain the chemical network 
that leads to water formation. A significant step in the under- 
standing of water formation was achieved thanks to the detec- 
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tion with APEX of HOOH in Oph A ( iBergman et"al~ll2011bl) . 
This detection has been carefully modeled with an astrochem- 
ical model combining reactions in the gas phase and on the 
grain sur faces based on the HME (Hybrid Moment Equation) 
method (iDu & Parisd 1201 ll). applied to the conditions of the 
Oph A cloud dDu et al.ll2012l) . The model successfully repro- 
duces simultaneously the abundances of the observed gas-phase 
HOOH, 2 , H 2 CO and CH3OH. The model also predicts that the 
hydroperoxyl radical (HO2, or O2H), the direct precursor lead- 
ing to HOOH via hydrogenation, has a significant abundance, 
which makes its detection with current radio telescopes possible. 
Hydroperoxyl has not yet - to our knowledge - been detected in 
the interstellar medium in either its solid or its gaseous form. 

We aim here at detecting this radical towa rd the SMI core 
of the p Oph A cloud, where HOOH and 2 dLiseau et al .1120121) 
were detected. The paper is organized as follows. The obser- 
vations are presented in Sect. 2, and the results in Sect. 3. We 
discuss the implications of this detection in Section 4. 

2. Observations 

The HO2 radical is a light asymmetric rotor with an unpaired 
electronic spin. This causes each asymmetric rotor level (except 
Oo,o) to be split into a doublet of states. This is illustrated in 
the HO2 energy level diagram in Fig. Q] where all states below 
50 K are shown. It is based on the JPL catalog dPickett et alj 
1998 ), which us e s mm and submm spectr oscopy data from 
Beers & Howard! (1 19751) . ISaitol (1 19771) . and ICharo & de Lucial 
(1982]). Also, the hyperfine splitting is large enough to be de- 
tectable. The spectrum of HO2 exhibits both a- and b-type tran- 
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sitions with associated dipole mo ments [i a = 1.41 D and ^ = 
1 .54 D (ISaito & Matsumuralll980h . 

We targeted lines of the HO2 radical at 130.3 GHz and 
260.6 GH z (see Table [D, which were observe d in the laboratory 
bv lSaitol (fl977l) and lCharo & de Lucial d 1982b . respectively. We 
used the frequencies tabulated in the JPL database, which cor- 
respond to the measured values for these lines. The lo,i — > Oo,o 
lines around 65 GHz are poorly suited for a ground-based search 
because they are heavily affected by the strong atmospheric ab- 
sorption due to O2. The 3o ; 3 — > 2o,2 lines (at 195 GHz) fall in a 
frequency region where few receivers are currently in operation. 
They are also somewhat affected by the atmospheric absorption 
from the 183 GHz water line. Likewise, the 5o,5 — > 4o,4 lines 
are affected by the atmospheric 325 GHz water line. Hence, the 
targeted lines at 130 and 260 GHz are the best candidates for a 
detection in a low-excitation source (for which most of the pop- 
ulation is likely to reside in the lower K a — levels). 
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Fig. 1. Energy level diagram of HO2 below 50 K. Only the a- 
type transitions listed in the JPL catalog were included and are 
drawn as arrows. Red arrows indicate observed transitions. The 
fo-type transitions (connecting K a = and 1 states) were omitted 
for the sake of clarity. Next to each level the quantum number 
designation is listed. 



All observations targeted the SMI position at 
a(2000)=16 /, 26"'27 s .20 , (5(2000)=-24°24'04". The IRAM 
30m observations were carried out on Dec. 5, 2011, under 
average weather conditions with a precipitable water vapour 
(PWV) column of 5 mm. The EMIR receiver dCarter et al.l 
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Fig. 2. Lines observed with the IRAM telescope. The red bars 
show the expected fre quencies for the different lines, as mea- 
sured in the laboratory (Saito fT977|). The blue bars s how the cal- 
culated frequencies from Charo & de Lucial (Il982l) . The black 
vertical line in the lower panel shows the position of the unde- 
tected F = 1 — > line. 
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Fig. 3. Lines observed with the APEX telescope. The data were 
smoothed to a resolution of 0.18 kms~'. The red bars show 
the expected f requencies for the differe nt lines, as measured in 
the laboratory (ICharo & de Lu cia 1982). The blue bars show the 
computed hyperfine frequencies (see Table[T}. 



|20H was tuned in both polarizations to 130.3 GHz. The 
angular resolution of the observations is 19" at this frequency. 
The receiver was connected to FTS spectrometers, which were 
used in their highest resolution mode, providing a 50 kHz 
resolution. The focus was optimized on Saturn at the beginning 
of the observations. The pointing was checked locally every 
hour, either on 1741-038, 1622-297, or 1730-130, depending on 
the time of the pointing, and was found to be better than 3'.'8 
(rms). For the conversion to the T m \, scale, we adopt a main 
beam efficiency of 0.76, as extrapolated from the values given 
on the IRAM website^. 

The APEX observations were carried out on Oct 10 and 16, 
2011, under variable weather conditions (PWV =2.3 - 3.7 mm 
on Oc t 10, PWV = 0.6 m m on Oct. 16). The APEX1 re- 
ceiver (Vassilev et al. 2008) was tuned to 260.25 GHz and con- 
nected to the XFFTS backend, leading to a velocity resolution 
of 0.09 km s _1 . The angular resolution of the observations at this 
frequency is 24". Because no planet was available during the 



http://www.iram.es/IRAMES/mainWiki/Iram30mEfficiencies 
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Table 1. The HO2 and SO2 lines detected in the same setups. 



Line 


Frequency 
(GHz) 


Quantum numbers 
Rotation /„ — > J t 


f „ -» F, 


(s" 1 ) 




E u 
(K) 


Su 


J T mb dv 
(mKs" 1 ) 


Av 
(km s- 1 ) 


(km s- 1 ) 


1 

2 
3 


130 26007 
130.25813 
130.46741 
130.46368 


2o 2 

2(1.2 
2o,2 
2o,2 


-> lot 
-» lo,i 
-> lo.i 
-» lo,i 


5/2^3/2 
5/2->3/2 
3/2->l/2 
3/2->l/2 


3— >2 
2-»l 
2-»l 
l-»0 


2.07x10" 
1.85x10" 
1.76x10" 
1.16x10" 


-5 
-5 
-5 

-5 


9.4 
9.4 
9.4 
9.4 


7 
5 
5 
3 


74+7 
54±8 
45±4 
< 26 


66+0 07 
0.65±0.11 
0.54±0.06 


4.6+0.1 
4.6±0.1 
3.0+0.1 


4a 
4b 
5a 
5b 


260.56614 
260.56566 
260.77030 
260.76967 


4 ,4 
4().4 
4o,4 
4 ,4 


-» 3o,3 
-» 3 ,3 
-> 3 ,3 
-> 3 ,3 


9/2^7/2 
9/2^7/2 
7/2^5/2 
7/2->5/2 


5^4 
4^3 
4^3 
3->2 


1.86x10" 
1.81x10" 
1.80x10" 
1.72x10" 


-4 
-4 
-4 
-4 


31.3 
31.3 
31.3 
31.3 


11 

9 
9 
7 


} 115+14 
} 61+13 


} 0.88±0.13 
} 0.91±0.24 




SO, 

so 2 


129.51481 
131.01486 


10 2 ,8 

12 W i 


-» 10l,9 
-> 12 .12 














199+6 
168±5 


0.64±0.02 
0.68±0.03 


3.8+0.1 
3.7+0.1 



Notes. For lines 1 to 3 and SO2, the frequencies are those measured in the laboratory and tabulated in the JPL database. The vi a of the lines are 
therefore computed based on these frequencies. The hyperfine structure frequencies of lines 4 and 5 were computed from the energies tabulated 
in the JPL database (http://spec.jpl.nasa.gov/ftp/pub/catalog/archive/c033001.egy I. Uncertainties on fluxes, linewidths and vj sr are lcr. The upper 
limit on the F = 1 — > line is the flux of the bump at the expected frequency (see Fig.0. 



observations, the focus was optimized on IRAS 15 194-51 15 by 
maximizing the CO(2-l) line flux. The pointing was checked ev- 
ery hour by observing CO(2-l) on RAFGL1922, and was found 
to be always better than 3". 



3. Results 

Fig.|2]displays the three lines targeted with the IRAM telescope, 
and Fig. [3] shows the two lines targeted with APEX. Over the 
wide bandwidth observed with the FTS, no obvious other line is 
present, except for two SO2 lines, which allow us to confirm the 
vi sl (3.8 kms" 1 ) and line width of the expected lines (see Table 

ID- 

We detected three HO2 lines with the IRAM telescope with a 
high signal-to-noise ratio, whose linewidths are consistent with 
the linewidth measured for the SO2 lines, but whose frequen- 
cies do not co i ncide exactly with those measured in the labora- 
tory by ISaitol d!977l) . shown as red vertic al bars in Fig. |2] (as- 
suming a vi sr of 3.8 kms" 1 ). ISaitol d 19771) estimated an uncer- 
tainty on the laboratory frequency of .2 MHz (i.e. 0.5 kms" 1 
at 130 GHz). ICharo & de Luc ia (1982) measured the frequency 
of other lines between 150 and 550 GHz, and used the full 
dataset available at that time (including Saito's lines) to fit the 
HO2 Hamiltonian. They found in particular a significan t dis- 
crepan cy (0.37 MHz) between the frequency measured by [Saitol 
dl977l) and the calculated frequency from the Hamiltonian for 
the two lines at 130.26 GHz, which they excluded from their 
fit. The blue bars i n Fig . [2] show the calculated frequency from 
|Ch aro & de Lucia ( 1982). The first two detected lines coincide 
very closely with their theoretical frequencies. When adopting 
these new frequency values, the vi sr of the lines become 3.8 and 
3.7 km s" 1 , respectively, which is fully consistent with the source 
velocity as determi ned by SO2 and sev eral other molecules in 
the SMI source (Bergman e t al.ll201 lah . The third line is still 
slightly offset, with a vi sr of 3.4 km s~' . 

Because the source is relatively line-poor, we can assign the 
three lines to HO2 with a high confidence level. An additional 
line (at 130.46368 GHz) is only tentatively detected, and we 
consider in the following only an upper limit on its flux. No 
blend with other plausible molecules was found at these frequen- 



cies in either the C DMS dMuller et al.ll200ll) or the JPL catalog 
dPickettet al.lll998l) . 

We detected two additional lines of HO2 with the APEX 
telescope. Both lines have an unresolved hyperfine structure. A 
Gaussian fit of these two lines leads to a broader line width than 
for the other lines (see Table [TJ, but this may be because of the 
hyperfine splitting, which for these lines corresponds to about 
0.6 -0.7 km s" 1 . Here again, no blend with any other plausible 
molecule is found. 

In summary, we report here the detection of five different 
spectral features assigned to seven different lines of HO2. 

4. Analysis and discussion 

We use the rotation diagram method to derive the abundance of 
HO2, which assumes on the assumptions that the population dis- 
tribution can be described by a single temperature, and that lines 
are optically thin. The rotation diagram obtained with the five 
detected spectral features is presented in Fig. |4] We assume for 
lines 4 and 5 that the measured flux is coming from the two un- 
resolved hyperfine components proportionally to A a ; x g u , as ex- 
pected for two optically thin lines originating from two levels 
of the same energ y. We assume a source size of 24", as for hy- 
drogen peroxide (Bergman et al.ll201 lbl) . and as derived by the 
analysis of Be rgman et al.l (1201 lab . We derive a rotational tem- 
perature of 16 + 3K, slightly lower than the rotational tem per- 
ature derived for HOOH (22 + 3 K, IBergman et ai1l201 lbl) . but 
marginally consistent with it within the error bars. The derived 
HO2 column density is (2.8 + 1.0) xlO 12 cm" 2 . This column 
density results in optically thin lines. Forcing the rotational tem- 
perature to be 22 K as for HOOH, we derive a column density 
for HO2 of 3.0 x 10 12 cm" 2 , very close to the column density 
derived with the fitted T mt of 16.4 K. 

The relatively high A u \ of the 260 GHz lines (and the even 
higher spontaneous rates of the b-type transitions) indicate that 
the excitation could b e somewhat subtherma l for the density in 
SMI (10 5 - 1 6 cm 3 . IBergman et al.ll20 1 1 ah . This could imply 
that our derived column density is too high. To estimate this, we 
simply truncate the partition function such that the summation 
only takes place over K a — states below 40 K and compare the 
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Fig. 4. HO2 rotational diagram. A source size of 24" is assumed. 
The triangles correspond to detected lines. The square corre- 
sponds to the upper limit of the flux on the undetected line at 
130.464 GHz. 



result to the full partition function. We then find that the trun- 
cated part holds about 70% of the level populations at T — 16 K. 
Therefore, in the unlikely case of extreme subthermal excitation 
(no population in K a > states and in K a = states above 
4o,4), we may explain the observations by adopting a 30% lower 
HO2 column density. This lower column density falls within the 
errors and we conclude that subthermal excitation will have a 
minor impact on our derived column density. 

From the H7CO a nd CH3OH analysis of the SMI core, 
Berg man et al.l (1201 lab determined an H2 column density of 
3 X 10 22 cm" 2 . The abundance of H0 2 i s therefor e ~10~ 10 . 

The gas-grain chemical model of iDu et al.l d2012l) . based 
on the detection of HOOH and other molecules, predicts that 
the HO2 / HOOH abundance ratio in p Oph A should be about 
three. In this model, at T ~ 20 K, HO2 is mainly formed on dust 
grains through the barrierless reaction O + OH at an early stage 
(< 5 x 10 4 yr), and through the reaction H + O2 at a later stage. 
Here, we find observationally H02/HOOH~ 1, which is con- 
sistent with the model, taking into account the excitation uncer- 
tainty. This result is therefore an additional validation of our un- 
derstanding of water formation on dust grains. 

Motivated by our detection, we inspected previous spec- 
tral surveys for unidentified lines at the frequencies of 
the HO ? molecule . No line is detected in the TIMASSS 
survey (ICaux et al. 



2011b toward the l ow-mass protostar 



Nummeli n et al.l (Q~998) reported (in their 



IRAS 16293-2422. 
Table 46) two unidentified lines in their Sgr B2 survey that 
might correspond to HO2 frequencies, U260568 and U260771 
toward Sgr B2(N) and U260569 toward Sgr B2(M). Because 
the lines have typical widths of 13kms~', line blending is a se- 
vere limiting factor for a secure line assignment. The Sgr B2 
IRAM 30 m line survey coverage does not include the lines at 
130 GHz (Belloche, priv. comm.). A 2mm unpublished sur- 
vey from B. Turner made with the NRAO 12m telescope be- 
tween 1993 and 1995 covers the two 130.2 6 GHz lines (but not 
the 130.46 GHz lines, Rem iian etalJ l2008). This doublet unfor- 
tunately falls very close to the 7=3-2 SiO emission/absorption 
feature at 130268.7 MHz, making it very difficult to discern 
any weaker HO2 lines. Assuming that the unid entified lines of 
the Sgr B2(N) survey by Nummel in et al.l (1 19981) are a ssigned to 
H0 2 , and assuming T ro i — 50 K (Nummelin et al. 2000), we find 
N(H02) < 1.1 x 10 14 cirT 2 , averaged in the 20" beam. Using the 
H2 column density derived by Nummelin et al. (3 X 10 24 cnT 2 ), 
this results in an upper limit on the fractional abundance of HO2 



of 4 x 10 11 , somewhat below the detected abundance toward 
p Oph A. 

The result of this work is one of the rare examples of the 
detection of a new molecule after it was predicted by chemi- 
cal modeling dDu et al.ll2012l). Other examples include the de- 
tection of CF + dNeufeld et alj ^6ey. c frf6tfvated also by (gas- 
phase) model predictions. The novelty here is that the detec- 
tion of HO2 was predicted by models describing the chem- 
istry on dust grain surfaces. This shows the increasing pre- 
dictive power of these models, thanks to the numerical re- 
finements develope d to correctly treat the stochasticity o f sur- 
face chemistry (e.g. lGarrod et al.ll2009l IDu & Parise l l2011l), and 
the r e cent numerous laboratory experiments (e.g. Oba et al.l 
120091: lloppolo et al.l l2010t ICuppen et al.l 120101) . The era when 
dust chemistr y could be describ ed as "the last refuge of the 
scoundrels" (Char nlev et al.l fl992) seems to be over. 
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